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摘 要 : 大 通 河流 域 地 处 青藏 高 原 东北 部 边缘 ,生态 环境 敏感 脆弱 ,开展 变化 环境 下 水 资源 演变 A 
因 研 究 对 保护 区 域 水 生态 环境 ,保障 水 生态 文明 建设 上 共有 十 分 重要 的 意义 。 采 用 线性 倾向 估计 、 
集中 度 、 集 中 期 有 序 聚 类 检验 、 小 波 分 析 等 统计 方法 ,分 析 了 流域 径流 的 年 际 变化 、 年 内 分 配 、 周 
期 和 突变 变化 特征 ,基于 累积 量 和 斜率 变化 率 法 和 双 累 积 曲 线 定量 评估 了 气候 因素 各 人 类 活动 对 
径流 变化 的 影响 。 结 果 表 明 :(1) 近 60a 大 通 河流 域 气候 暖 湿 化 明显 ,年 平均 气温 、 降 水 量 、 潜 在 
蒸发 量 增幅 分 别 为 0.42 C (10a)! F 8.9 mm: (10a) ', 5.6 mm- (10a)! ,年 径流 呈 减 少 趋势 ,倾向 率 
0.67X10° m°- (10a)' (2) 径流 集中 度 和 不 均 久 系数 呈 微 弱 下 降 趋 势 , 枯 季 径 流 增 加 趋势 明显 ,年 内 
分 配 更 趋 于 均匀 ,集中 期 有 推迟 趋势 , 延 诬 速率 为 3.0d':(10aj (3) 年 径流 在 44 a 左右 尺度 上 周期 
震荡 明显 ,突变 发 生 在 1990 年 ,突变 后 径流 量 减少 3.52x10* m? ,流域 冰川 分 布 呈 减 小 趋势 ,植被 履 
盖 无 显著 变化 ,(4) 气候 、 人 类 活动 对 大 通 河 径流 减少 的 贡献 率 分 别 为 -17.7% 和 和 117.7%, 降 水 量 是 
流域 来 水 的 主要 补给 来 源 , 跨 流 域 调 水 是 引起 径流 减 小 的 主要 驱动 因素 。 
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大 通 河 地 处 青藏 高 原 东北 边缘 ,生态 环境 脆 
弱 , 是 我 国 西北 地 区 重要 的 生态 过 渡 带 , 随 着 水 能 
水 资源 开发 等 人 类 活动 的 加 剧 ,河流 水 文 过 程 、 水 
生态 环境 受到 严重 影响 。 另 外 , 受 全 球 气 候 暖 湿 化 
影响 ,流域 降水 径流 等 关键 水 循环 要 素 也 发 生 了 
显著 变化 ,人 类 活动 加 剧 的 水 循环 不 仅 会 影响 到 水 
资源 生态 环境 ,还 将 使 下 游 地 区 的 水 资源 安全 面临 
诸多 风险 ,因此 开展 归 因 分 析 , 研 究 其 水 文 变 化 机 
制 对 流域 水 资源 科学 管理 和 保护 利用 具有 重要 意 
义 。 当 前 ,统计 分 析 法 .水 文 模型 法 和 弹性 系数 法 
等 是 定量 识别 径流 变化 归 因 的 主要 方法 "”。 吴 恒 
卿 等 ”采用 分 布 式 水 文 模型 SWAT 对 大 通 河流 域 径 
流 序列 进行 了 模拟 。 白 雁 鱼 等 ”基于 Budyko 假设 
的 气候 弹性 系数 法 分 解 了 降水 潜在 蒸发 和 其 他 因 
素 对 径流 变化 的 影响 分 量 。 王 大 超 ” 采用 双 累 积 曲 
线 得 出 人 类 活动 在 1995 一 2007 年 .2008 一 2013 年 对 
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大 通 河 下 游 径流 的 影响 紊 分别 为 92.5% 和 135.4%。 
水 文 模型 参数 率 定 和 验证 过 程 复 杂 ANE PET, 
基于 Budyko 假设 的 弹性 系数 法 主要 在 年 尺度 上 定 
量 区 分 影响 因素 对 径流 变化 的 影响 ,部 分 学 者 在 归 
识别 时 因 方 法 本 身 .参数 选取 、 计 算 时 段 不 同 往 
往 得 出 不 同 的 结论 。 以 往 研 究 倾向 于 采用 一 种 方 
法 来 估计 径流 变化 的 贡献 率 ,结果 缺少 不 同方 法 的 
对 比 验 证 , 且 对 径流 影响 因素 分 析 多 为 定性 描述 。 
王 随 继 等 "提出 的 累积 量 斜 率 变化 率 法 所 需 数据 简 
单 ,累积 量 在 一 定 程度 上 消除 了 人 为 干扰 以 及 实测 
数据 年 际 波动 的 影响 ,被 诸多 学 者 应 用 于 径流 归 
研究 中 。 笔 者 基于 最 新 的 水 文 气象 资料 ,首次 采用 
累积 量 斜率 变化 率 法 对 大 通 河 流域 径流 变化 的 气 
候 变 化 、 人 类 活动 等 影响 因素 进行 定量 评估 * ,同时 
利用 双 累 积 昌 线 法 对 结果 进行 分 析 验 证 ,以 揭示 流 
域 水 文 过 程 对 环境 变化 的 响应 关系 ,人 研究 方法 和 成 
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果 可 为 流域 综合 开发 治理 提供 科学 依据 。 
1 研究 区 概况 


大 通 河 为 涅 水 最 大 的 一 级 支流 ,发 源 于 青海 省 
天 峻 县 托 勒 南山 ,地 理 位 置 介 于 98°30'~103°15'E、 
36°30'~38°25'N 之 间 ,河源 海拔 4520 m ,干流 全 长 
572 km, ,河口 海拔 1727 m, 河 道 平均 比 降 4.19%。 河 
WERKIE LYE, IK 308.4 km, KARA, Œ 
树枝 状 ,主要 支流 有 葛 日 曲 、 江 仓 曲 、 永 安河 等 7 
条 ,主要 分 布 在 右岸 ; 角 大 滩 至 连城 为 中 游 段 , 河 
长 223.4 km ,水 系 呈 羽 状 ,主要 支流 有 白水 河 .老虎 
沟 、 讨 拉 沟 等 ;连城 以 下 为 下 游 段 , 河 长 40.2 km 
(图 1)。 大 通 河 流域 属 内 陆 高 寒气 候 区 ,气候 特点 
为 冬 长 夏 短 冰冻 期 长 .气候 寒冷 .年 降水 较 多 、 雨 
热 同 季 , 流 域内 年 平均 气温 -0.3~8.0 ,极端 最 高 气 
温 35.8 ,极端 最 低 气温 -34.1 °C, BAPE IY Bek E 
500.6 mm, 


2 数据 与 方法 


2.1 数据 来 源 

大 通 河 出 口 处 布设 有 享 堂 水 文 站 ,控制 流域 面 
积 15126 km’, HEA H O 1.9 km ,该 站 设立 于 1939 年 
10 月 ,1940 年 开始 观测 ,经 黄 委 整编 后 的 完整 连续 


资料 自 1950 年 至 今 。 流 域 和 周围 10 处 气象 站 的 逐 
日 降水 量 .气温 、 相 对 湿度 日照 时 数 、 风 速 等 资料 
来 源 于 中 国 气象 数据 网 (http:/data.cma.cn ) , 面 降 水 
量 PERRE .平均 气温 采用 泰 森 多 边 形 法 计算 ， 
考虑 各 站 点 资料 的 时 间 序 列 和 完整 性 ,水 文 气象 要 
素 分 析 统 一 采用 1961 一 2020 年 系列 。 冰 川 分 布 
和 地 表 植 被 指数 产品 "来 源 于 国家 青藏 高 原 科学 
数据 中 心 (http://data.tpdc.ac.cn) ,冰川 分 布 产 品 原始 
基础 数据 为 1978 一 2019 年 祁连山 全 境 的 Landsat 系 
列 影像 ,时间 分 辩 率 为 Sa, 空 间 分 辩 率 为 30 m, 用 于 
冰川 变化 对 流域 径流 影响 定量 估计 。 地 表 植 被 指 
数 产 品 空间 分 辩 率 为 30 m, 采 用 最 大 值 合成 方法 实 
现 对 地 表 月 度 归 一 化 植被 指数 (Normalized differ- 
ence vegetation index,NDVJT) 产 品 的 合成 ,进而 计算 
植被 覆盖 度 (Fractional vegetation cover, FVC). 
2.2 研究 方法 
2.2.1 HAZE KH Penman-Monteith 公式 上 计 
算 流 域 潜在 蒸发 量 : 
0.408A(R,- G) + py a ~e,) m 
A+y( +0.34u,) 
式 中 :ET 为 潜在 蒸发 量 (mm*d"'); A 为 饱和 水 汽 不 
曲线 斜率 ; R, 为 地 面 净 辐 册 (MJ"m …d ); 6 为 土壤 
PU at (MJ-m?+d'); y AMER TT ae Be (kPa C); 
7 为 日 平均 气温 (C); w, 为 2m 高 处 风速 (ms '); 


ET= 


图 例 
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图 1 大 通 河流 域 水 系 及 站 点 分 布 


Fig. 1 Distribution of water System and stations in the Datong River Basin 
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e. 为 饱 合 水 汽 压 (kPa); e, 为 实际 水 汽 压 (kPa)。 
2.2.2 趋势 分 析 及 突变 检验 采用 线性 倾向 估计 于 
分 析 气 象 水 文 要 素 的 变化 趋势 。 根 据 变异 诊断 分 
析 方 法 ,采用 过 程 线 、 累 积 距 平 进行 年 径流 量 的 初 
步 诊 断 , 确 定 序列 是 否 发 生变 异 , 通 过 Mann-Kend- 
all AP ERE 、Pettitt 检 验 "详细 诊断 识别 径流 的 
突变 年 份 。 采 用 集中 度 、 集 中 期 .不 均匀 系数 “等 
指标 反映 河川 径流 年 内 分 配 均 匀 程 度 。 

(1) 线性 倾 回 估计 

H y 表示 样本 为 n 的 变量 , x 为 y, 所 对 应 的 时 
间 ,建立 y; 与 x 的 线性 回归 方程 : 

y,=ax,+b, i1=1,2,.…,n (2) 

式 中 ;a 为 回归 系数 ;。 为 回归 常数 ;a 、b 值 可 采用 
最 小 二 乘法 估计 。 


式 中 ;7 为 径流 量 序列 的 平均 值 ; x 为 时 间 序 列 序 
写 的 平均 值 。 可 求 出 时 间 x 与 变量 y, 之 间 的 相关 
系数 7 


(2) 有 序 聚 类 检验 
其 实质 在 于 求 算 最 优 分 割 点 , 即 寻 找 序列 突变 
前 后 系列 离 差 平方 和 综 和 最 小 的 点 。 


v= Pea), = (5) 


t=T+1 


SUV, V, 分 别 为 分 割 点 前 后 离 差 平方 和 ;7 为 
可 能 分 割 点 ; x, 为 序列 值 ; 元 为 突变 点 7 前 的 序列 
均值 ; 元 ,为 突变 点 7 后 的 序列 均值 。 

S,n=V,+V,_, (6) 

式 中 : S(t) 为 总 离 差 平方 和 ,其 取 最 小 值 时 对 应 的 
T 即 为 最 优 分 割 点 。 
(3) Pettitt 检 验 

Pettitt 检验 为 非 参数 检验 法 ,物理 意义 清晰 ,能 

较 好 地 识别 水 文 气象 时 间 序 列 中 的 突变 点 。 对 于 具 

有 样本 量 为 n 的 时 间 序 列 x ,可 以 构造 一 秩序 列 为 : 

oye k=2,3,.…,n (7) 


其 中 : 


十 1 Xi > x; 
Se ds (8) 
1,0 < 
式 中 : se 为 统计 量 ; k 为 样本 总 数 ;P 为 显著 性 水 
平 ; x, x, 为 相同 时 间 序 列 的 2 个 样本 值 ; r 为 符号 
消 数 , 即 第 i 数值 大 于 j 时 刻 数值 个 数 的 累计 数 。 
Pettitt 法 是 直接 利用 秩序 列 来 检验 变异 点 的 ， 
Fe t WAN Æ : 
k, =max|s, p|, k=2,3,.…,n (9) 
式 中 :为 最 大 统计 量 ,对 应 的 t 点 为 变异 点 ,计算 
统计 量 P=2exp|-6k,*(n’ +n’)] “4 P<0.5 HY, WIG 
检验 出 的 变异 点 在 统计 意义 上 是 显著 的 。 
2.2.3 小 波 周期 分 析 小 波 分 析 "! 是 在 傅 里 叶 变 化 
基础 上 发 展 起 来 的 一 个 数学 分 文 ,能 清晰 揭示 出 隐 
含 在 时 间 序 列 中 的 多 种 变化 周期 ,被 广泛 应 用 于 水 
文 气象 序列 周期 识别 。 小 波 方差 可 识别 各 种 尺度 
的 震荡 强 弱 和 周期 变化 特征 ,方差 值 越 大 ,表明 对 
应 尺度 的 周期 越 明 显 。 对 于 小 波 函 数 VD ,径流 时 
FR x(t) L (R) 的 连续 小 波 变换 为 : 
w ,(a,b)=— xp += Jae 10 
(ab) eles (10) 
式 中 : wr 为 小 波 变换 系数 ; a 为 尺度 因子 ;5 为 
平移 因子 ; R 为 实数 域 ; x) 为 径流 序列 ; Wx) 为 
2.2.4 径流 变化 归 因 识别 
(1) 累积 量 斜 率 变化 率 法 
累积 量 冬 率 变化 率 法 是 在 基准 斯 和 突变 期 分 
别 建立 影响 因子 与 年 份 之 间 的 线性 方程 ,从 而 通过 
和 斜率 变化 率 的 比值 计算 出 影响 因子 的 贡献 率 。 累 
积 降水 量 .潜在 蒸发 量 和 径流 量 的 斜率 变化 率 分 
别 为 : 


Sp = (Ky, ~ K,,)/Kp, ( 11 ) 
Ser = (Kin ~ Kier)/Ker, ( 12) 
SR =( 人 mu 一 大 人 人。 (13) 


中 : Sp. Sup. Sy 分 别 为 累积 降水 量 EEA At. 
流量 变化 期 相对 于 天 然 时 期 的 变化 率 (%); K,、 

发 量 径流 量 突变 前 后 线性 关系 式 的 斜率 mm'a-)， 

正 值 ( 负 值 ) 表 示 和 斜率 增 大 ( 减 小 )。 

通常 情况 下 径流 量 与 降水 量 呈 正 相关 ,潜在 蒸 
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C, =S,/8, X 100% (14) 
Cir =—S wr! Sy X 100% (15) 
式 中 : Cy. Ca 分 别 为 降水 量 EE ZS ENE 
化 的 贡献 率 。 
人 类 活动 对 径流 变化 的 贡献 率 Cu : 
Cy=1-C,— Cr 
(2) 双 累 积 曲 线 法 
采用 双 累 积 曲 线 法 对 大 通 河流 域 径 流 变化 归 
因 结 果 进 行 验证 。 双 累积 曲线 "是 水 文学 中 销 用 
来 反映 降水 径流 系统 变化 的 特征 曲线 ,天然 情况 下 
累积 降水 量 和 累积 径流 量 应 为 一 条 斜率 固定 的 直 
线 , 当 斜率 发 生 显著 变化 时 说 明 人 类 活动 对 径流 的 
影响 持续 增加 ,该 年 份 可 判定 为 突变 点 。 采 用 基准 
期 降水 径流 关系 和 变化 期 面 降水 量 计算 仅 受 降水 
影响 的 理论 径流 量 ,与 基准 期 的 径流 量 差 值 为 降水 
变化 的 影响 量 ,扣除 降水 因素 导致 的 径流 影响 量 即 
为 人 类 活动 对 径流 变化 的 影响 量 。 


(16) 
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3 结果 与 分 析 


年 际 变化 趋势 分 析 
近 60 a 大通 河流 域 气候 呈 暧 湿 化 方向 发 展 ( 表 1、 
图 2a~b) ,年 平均 气温 表现 出 非常 明显 的 上 升 趋势 ， 
增幅 0.42 C- (10a) 明显 高 于 全 球 [0.26 C- (10a) ]™, 
通过 0.001 的 显著 性 检验 ;年 降水 量 呈 增加 趋势 , 倾 
向 率 为 8.9 mm'.(10aj…。 淤 在 蒸发 量 呈 增加 趋势 , 增 
幅 5.6 mm (10a)'。 年 径流 量 呈 减少 趋势 ,变化 速率 
0.67x10” m…(10a) ,最 大 值 50.19x10* m (1989 年 )， 
最 小 值 19.95x10: m*(2015 年 ), 极 值 比 2.52, 变 差 系 
数 为 0.20, 年 降水 量 潜在 蒸发 量 、 年 径流 量 均 通过 
0.05 的 显著 性 检验 。 
3.2 ”径流 年 内 变化 特征 

大 通 河流 域 径流 年 内 分 配 的 集中 度 、 不 均匀 系 
数 (图 3a~b) 呈 微弱 下 降 趋势 ,径流 年 内 分 配 趋 于 均 
匀 , 汛 期 (5 一 10 月 ) 径 流量 占 全 年 比重 呈 减 小 趋势 ， 
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表 1 大 通 河流 域 各 要 素 趋势 检验 


Tab.1 Trend test of elements in the Datong River Basin 


水 文 气象 变化 ie a Kendall a - Spearman asa 
EA TEK EA TEK EA IEK 
BR ae 4 (a=0.05 š tan=2.0) ( a=0.05 $ Un=1.96 ) (a=0.05 š tan=2.0 ) 
年 平均 气温 增加 16.30 趋势 显著 7.74 趋势 显著 16.70 趋势 显著 
年 降水 量 增加 2.58 趋势 显著 299 趋势 显著 2.25 趋势 显著 
潜在 蒸发 量 增加 1.83 趋势 不 显著 1.08 趋势 不 显著 1.14 趋势 不 显著 
年 径流 量 减少 172 趋势 不 显著 1.53 趋势 不 显著 1.51 趋势 不 显著 
注 :7T.0 为 统计 量 。 
一 一 一 年 径流 量 ------ Salil 一 一 一 潜在 蒸发 量 ------ 5a 滑动 
一 一 年 降水 量 ------ 趋势 线 一 一 年 平均 气温 ------ 趋势 线 
80- (a) 年 降水 量 与 径流 量 700 8r (b) 年 平均 气温 与 潜在 蒸发 量 1100 
70 600 í y=0.5612x—192.64 inii 
Pal 6 R=0.0543 
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图 2 各 要 素 序列 变化 趋势 


Fig. 2 Sequence change trend of each element 


2 期 王 淑 芝 等 : 


而 枯 季 (11 月 至 翌年 4 月 ) 径 流量 显著 增加 ,反映 出 流 
域 生 态 环境 好 转 , 蕾 水 和 水 源 涵养 能 力 逐 步 增强 。 随 
时 间 推 移 ,集中 期 有 推迟 趋势 ,延迟 速率 为 3.0 de 
(10aj” ,2011 一 2020 年 最 大 径流 量 出 现时 间 较 20 世 纪 
60 年 代 平均 延迟 18 d。 从 季节 分 析 , 夏 季 径 流量 ( 占 
比 48.1% ) 远 大 于 其 他 季节 ,呈现 夏季 (6 一 8 月 ) > 秋季 
(9 一 11 月 )> 春 季 (3 一 5 月 )> 冬 季 (12 至 翌年 2 月 )。 
3.3 径流 周期 性 分 析 

享 笛 站 小 波 系 数 实 部 等 值 线 和 小 波 方法 曲线 
如 图 4a~b 所 示 ,大 通 河 年 径流 主要 存在 12 a、24 ay 
44 a 的 周期 震荡 。3 个 时 间 尺 度 中 44 a 峰值 最 明显 ， 
说 明 周 期 震荡 最 强 , 为 第 一 主 周期 ,其 次 为 24a 和 
12 a。 在 44 a 左右 的 特征 时 间 尺 度 下 ,年 径流 经 历 
了 2 次 枯 - 丰 转化 ,当前 的 丰 水 期 直至 2021 年 等 值 
线 仍 示 闭合 ,表明 偏 丰 状态 仍 将 持续 。 
3.4 径流 突变 特征 分 析 

从 大 通 河 流域 降水 量 、 径 流量 累积 距 平 曲线 
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207 


(图 $) 初 步 判断 ,年 径流 量 在 1990 年 左右 发 生 突变 ， 
由 丰 水 期 转 为 枯水期 。Mann-Kendall 检验 表明 突变 
发 生 在 1995 年 ;有 序 聚 类 检验 曲线 离 差 平方 和 在 
1990、1994 年 达到 最 小 值 ,出 现 明显 的 突变 ;Pettitt 
检验 中 1990 年 径流 突变 点 统计 量 P=0.063<0.5 , 为 
有 效 突变 点 。 综 合 各 检测 方法 结果 ,确定 大 通 河 年 
径流 在 1990 年 发 生 突变 。20 世 纪 90 年 代 末 大 通 河 
流域 水 电 开 发 持续 加 快 ,1994 年 引 大 入 秦 调 水 工程 
建成 通 水 是 引起 流域 径流 突变 的 主要 原因 ,后 外 调 
水 量 逐 年 增加 ,流域 径流 量 明显 减少 。 相 较 基 准 
期 ,变化 期 径流 量 减 少 3.52x10: mi, 减 幅 12.1%; 降 
水 量 增 加 14.5 mm, 增幅 3.0% ; 潜在 蒸发 量 增加 
11.1 mm ,增幅 1.2%. 
35 径流 变化 定量 评估 

绘制 大 通 河流 域 年 径流 量 、 降 水 量 和 湾 在 蒸发 
量 的 累积 曲线 ,对 累积 曲线 的 基准 期 .变化 期 分 别 
建立 累积 量 与 年 份 的 线性 关系 , 拟 合 方程 (图 6) ,各 


一 一 不 均匀 系数 一 一 枯 季 径流 量 ----- 趋势 线 


2.6r (b) 不 均匀 系数 、 枯 季 径 流量 -10.0 

ae y=0.0311x-56.611 fi (80 

ie R=0.1529 i g 
= | a a epee 6.0 = 
g 14H | 
we | 4.0 X 
长 1.0 向 

0.6 NE ee 20 Ë 

y=-0.002x+4.8515 

R?=0.0802 
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图 3 享 得 站 径流 量 指标 年 际 变化 


Fig.3 Interannual variation of runoff indicators at Xiangtang Station 
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(b) 小 波 方差 曲线 
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图 4 径流 小 波 周期 分 析 


Fig. 4 Period analysis of runoff wavelet 
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(a) 降水 量 、 径 流量 累积 距 平 


径流 量 累积 距 平 /108m 
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注 :UF 为 时 间 序 列 统计 量 ;UB 为 逆序 统计 量 ;$, 为 总 离 差 平方 和 。 
图 5 大 通 河 流域 年 径流 突变 检验 
Fig. 5 Abrupt change of annual runoff in the Datong River Basin 


序列 在 1990 年 前 后 斜率 变化 率 有 明显 的 变化 特征 ， 
累积 降水 量 ,潜在 蒸发 量 的 斜率 变化 率 为 正 , 累积 径 
流 深 的 斜率 变化 率 为 负 ,说 明 降 水 量 BERRE 
列 呈 增加 趋势 ,径流 深 序 列 为 下 降 趋 势 。 由 式 (11)~ 
(13) 计 算 各 累积 序列 影响 期 的 斜率 变化 率 , 由 式 
(14)~(16) 计 算 气 修 变化 和 人 类 活动 对 径流 变化 的 
贡献 率 , 计 算 结 果 见 表 2。 结 果 表 明 , 气 候 变 化 中 降水 
量 增加 对 径流 量 减 少 呈 人 负责 献 ,贡献 率 -31.4%, 洪 在 
蒸发 量 增 加 对 径流 量 减 少 呈 正 贡献 ,贡献 率 为 
13.7% ; 人 类 活动 导致 径流 量 减 小 的 贡献 率 达 
117.7%。 因 此 人 类 活动 是 导致 流域 径流 变化 的 主导 
素 ,结论 与 杨 大 文 等 的 研究 结果 一 致 ,结果 可 信 。 

3.6 WHR HAA RIG iE 


i=ZA 
TA, 


理性 (图 7a~b) ,突变 后 的 降水 量 - 径 流量 关系 曲线 
和 斜率 明显 小 于 天 然 时 期 , 产 流 量 减 小 。 经 计算 ,大 
通 河流 域 1991—2020 年 径流 量 相 比 基 准 期 减 小 
3.52x10° mm ,气候 变化 导致 径流 量 增加 1.40108 mm， 
影响 率 为 -39.8% ,而 人 类 活动 导致 径流 量 减 小 
4.92108 m ,影响 率 达 139.8% , 即 人 类 活动 是 导致 
流域 径流 变化 的 主要 因素 ,与 累积 量 斜率 变化 率 法 
计算 的 贡献 率 接近 ,计算 结果 可 信 。 

3.7 径流 变化 影响 因素 分 析 

3.7.1 气象 因子 对 径流 的 影响 分 析 影响 径流 变化 
的 主要 气象 因子 包括 降水 .气温 和 蒸 散发 等 。 大 通 
河流 域 径 流 与 降水 的 年 内 分 配 规律 一 致 ,汛期 5 一 
10 月 径流 量 占 年 径流 量 的 80.4% ,降水 对 径流 表现 


B 
A; 


采用 双 累 积 曲线 验证 径流 变化 归 因 成 果 的 合 为 正 贡 献 ,是 流域 来 水 的 主要 补给 来 源 , 蒸 散发 则 
表 2 气候 变化 和 人 类 活动 对 大 通 河 流域 径流 影响 
Tab. 2 Impacts of climate change and human activities on runoff in the Datong River Basin 
径流 深 气候 变化 i 
+ — = 一 人 类 活动 
ale 斜率 变化 率 /% I r= 
ü AOOO O RR 变化 率 /% Cal o 斜率 IERO Cdo Cul% 
基准 期 (1961 一 1990 年 ) 187.2 z 483.4 一 一 921.8 = = = 
变化 期 (1991 一 2020 年 ) 168.1 -10.2 499.0 3.2 -31.4 934.7 1.4 13.7 117.7 


TE: Cp Cer, Cu A AREKE .潜在 蒸发 量 . 人 类 活动 对 径流 变化 的 贡献 率 。 
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图 6 大 通 河 流域 累积 年 径流 深 、 降 水 量 和 潜在 蒸发 量 与 年 份 的 拟 合 关系 


Fig.6 Fitting relationship between cumulative annual runoff depth, precipitation, evapotranspiration and 
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图 7 大 通 河流 域 降水 量 -径流 量 双 累 积 曲 线 及 其 关系 


Fig.7 Double-mass curve and relationship of precipitation-runoff in the Datong River Basin 


表现 为 负 贡 献 。 人 气温 升 高 一 方面 通过 影 
减少 河川 径流 量 , 男 一 方面 则 引起 流域 内 积 雪 和 冰 
川 消融 ,间接 增加 河川 径流 量 。 解 译 流域 内 冰川 面 
积 ( 表 3) 发 现 ,1986 一 2019 年 冰川 面积 整体 呈 减 小 


响 蒸 散发 ”趋势 ,1986 年 流域 冰川 面积 为 37.0 km’, 2005 年 减 小 
至 21.3 km’, 2015 年 以 后 面积 稳定 在 16.3$ km 左右 。 
绘制 年 径流 与 年 平均 气温 相关 关系 为 负 相 关 ,分 析 


流域 平均 气温 对 年 径流 表现 为 负 贡 献 。 因 此 ,气候 


表 3 大 通 河 流域 冰川 面积 及 植被 覆盖 度 统计 


Tab.3 Statistics of glacier area and fractional vegetation cover (FVC) in the Datong River Basin 


指标 1986 年 。 1990 年 ” 1995 年 ” 2000 年” 2005 年 ”2010 年 ”2015 年 ”2017 年 ”2018 年 。 2019 年 
冰川 面积 /km 37.0 29.5 25.0 223 21.3 29.5 16.3 16.3 16.4 17.3 
FVC 年 值 0.597 0.556 0.621 0.519 0.476 0.536 0.661 0.543 0.570 0.595 
FVC 生 长 季 均 值 0.723 0.744 0.762 0.688 0.688 0.699 0.827 0.727 0.727 0.726 
年 径流 次 /mm 191.3 196.3 171.3 158.4 189.5 143.0 131.9 193.0 175.3 198.8 
生长 季 径 流 深 /mm 114.7 97.6 73.0 71.1 83.2 73.9 69.4 64.3 56.9 72.5 


注 :FVC 为 植被 覆盖 度 。 
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条 件 中 降水 对 径流 变化 起 主导 作用 ,气温 SI 
冰川 面积 影响 径流 变化 但 不 显著 。 

3.7.2 植被 变化 对 径流 的 影响 分 析 大 通 河 流域 人 
烟 稀 少 , 开 发 利用 程度 低 ,土地 利用 类 型 以 林地 和 
草地 为 主 ,面积 达 70% , 短 时 期 内 可 认为 下 垫 面 地 
JÉ 地质 和 土壤 条 件 相 对 稳定 ,因此 下 垫 面 变化 更 
多 的 归 因 于 植被 功能 变化 。 植 被 变化 对 径流 的 作 
用 ,部 分 学 者 认为 植被 覆盖 度 (FVC) 增 加 会 减少 区 
域 径流 量 , 也 有 学 者 认为 FVC 增 加 会 提高 区 域 的 产 
流 能 力 。 从 表 3 可 知 ,1986 一 2019 年 流域 植被 变化 
趋势 不 明显 ,相关 关系 分 析 流 域 FVC 年 值 . 生 长 季 
均值 (6 一 8 月 ) 与 径流 相关 性 弱 。 考 虑 上 述 因 素 得 
出 1990 一 2020 年 大 通 河流 域 植 被 变化 对 径流 的 影 
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座 ) ,总 装机 容量 $66.97 MW ,分 布 在 中 下 游 河 段 , 建 
成 时 间 集 中 在 1990 一 2016 年 (图 8a) ,其 中 干流 23 
座 , 装 机 548.2 MW ,老虎 沟 等 8 条 支流 17 座 ,装机 
18.8 MW。 除 2014、2016 年 相继 建成 的 两 座 大 型 水 
库 纳 子 峡 和 石头 峡 水 库 对 流域 径流 洪水 有 调 蓄 作 
用 外 ,其余 径流 引水 式 电站 对 来 水 调 蓄 影响 较 小 ， 
水 库 运 行 对 河道 径流 的 影响 表现 为 对 天 然 来 水 进 
行 调 节 , 蕾 丰 补 村, 水 库 下 游 汛 期 径流 比重 、 集 中 度 
减 小 , 枯 季 径流 增加 (图 8b)。 

(2) 调 水 工程 

根据 大 通 河水 量 分 配方 案 , 大 通 河 设计 外 调 水 
量 14.43x10° mi’, 其 中 引 大 济 西 2.50x10 m’, 引 大 济 
1E 7.50x10° m, FI RAZR 4.43x10° ms BEA SLAF 


啊 较 小 。 西 一 期 引 硫 济 金 . 引 大 入 么 工程 相继 建成 发 挥 效 
3.7.3 水 利 工 程 对 径流 的 影响 分 析 益 , 流 域外 调 水 量 逐 年 增加 (图 9a)。 其 中 引 硫 济 金 


(1) 车 水 工程 
经 调查 ,大 通 河 流域 蓄 水 工程 主要 为 水 电站 ， 
无 汶 地 坝 工程 。 流 域内 已 建 水 电站 40 座 (含水 库 两 


工程 从 大 通 河 支流 硫磺 沟 调 水 至 甘肃 省 金昌 市 , 设 
计 引 水 规模 0.40x10 m, 2003 年 5 月 建成 通 水 ， 
2003—2021 年 平均 调 水 量 0.25x10 m°’, 2011 年 达到 
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图 8 水 电站 装机 容量 累积 及 不 同年 代 月 径流 占 比 


Accumulation of installed capacity of hydropower stations and distribution of monthly runoff in different years 
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图 9 调 水 工程 引水 量 及 享 符 站 天 然 年 径流 量 


Fig.9 Diversion of water diversion project and natural annual runoff hydrograph of Xiangtang Station 
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设计 调 水 量 。 引 大 入 秦 是 在 天 党 水 文 站 下 游 500 m 
将 大 通 河水 引 至 秦王 川 ,设计 引水 量 4.43x10 m, 
1994 年 10 月 工程 建成 通 水 ,1994 一 2021 年 平均 调 
水 量 1.91x10° m ,年 最 大 调 水 量 3.70x10s m(2021 
年 )。 引 大 济 涅 工程 目前 正 处 于 建设 阶段 ,2019 一 
2021 年 问 涅 水 上 游 宝 库 河 平均 调 水 1.17x10° m°, 
通过 对 享 堂 站 实测 径流 量 进行 还 原 , 经 肯 德 尔 秩 次 
检验 , 享 党 站 天 然 年 径流 呈 显 著 增加 趋势 (图 9b)。 
2019 、2020 年 流域 外 调 水 量 占 享 堂 站 径流 还 原 量 的 
90% 以 上 ,因此 跨 流域 调 水 是 导致 流域 径流 量 减少 
最 主要 的 关键 因素 。 


4 结论 


本 文 利用 大 通 河流 域 气象 水 文 等 数据 ,分 析 了 
近 60 a 径 流 演变 特征 ,运用 累积 量 斜 率 变化 率 法 和 
双 累 积 时 线 法 对 径流 变化 进行 归 因 分 析 , 得 出 主要 
结论 如 下 : 

(1) 大 通 河 流域 气候 暖 湿 化 明显 ,气温 、 降 水 
量 潜在 蒸发 量 增幅 分 别 为 0.42 °C (10a) 8.9 mm: 
(10aj #115.6 mm: (10a)"', 年 径流 量 呈 减少 趋势 , 倾 癌 
XH 0.67x10° m°- (10a), Æ 1990 FRERE , RA 
ALEID 3.52x10° m, IIR 12.19%. 

(2) 大 通 河流 域 径 流 集中 度 和 不 均 久 系数 呈 
微弱 下 降 趋 势 , 枯 季 径流 增加 趋势 明显 ,径流 年 内 
分 配 趋 于 均匀 ,集中 期 有 推迟 趋势 ,延迟 速率 为 
3.0 d- (10a)。 

(3) 大 通 河 年 径流 具有 12 a24 a N 44a EAH 
震荡 周期 ,在 44 a 左右 尺度 上 周期 震荡 最 为 显著 ,为 
第 一 主 周期 。 最 近 丰 水 期 的 小 波 等 值 线 还 未 闭合 ， 
表明 未 来 一 段 时 间 年 径流 将 持续 偏 丰 。 

(4) 气候 、 人 类 活动 对 大 通 河 径流 减少 的 贡献 
率 分 别 为 -17.7% 和 117.7% ,降水 是 流域 来 水 的 主要 
补给 来 源 , 跨 流域 调 水 是 引起 径流 减 小 的 主要 人 为 
因素 。 
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Attribution analysis of runoff change in the Datong River 
Basin, Qinghai-Tibet Plateau 
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2. Qinghai Province Flood and Drought Disaster Prevention Service Center, Xining 810001, Qinghai, China) 


Abstract: The Datong River Basin is located on the northeastern edge of the Qinghai-Tibet Plateau and is a sensi- 
tive and fragile ecological environment. It is of great significance to conduct research on the evolution and attri- 
bution of water resources in changing environments for the protection of the water ecological environment in the 
area and the construction of water ecocivilization. Statistical methods such as linear regression, concentration de- 
gree, concentration period, ordered clustering test, and wavelet analysis were used to analyze the characteristics 
of annual variation, seasonal distribution, periodicity, and abrupt changes in basin runoff. On the basis of the cumu- 
lative slope change rate method and double cumulative curve, the effects of climate factors and human activities on 
runoff changes were quantitatively evaluated. The results showed the following: (1) The climate in the Datong Riv- 
er Basin had warmed and humidified significantly in the past 60 years, with increases in average annual tempera- 
ture, precipitation, and potential evapotranspiration of 0.42 °C+(10a) ', 8.9 mm-(10a) ', and 5.6 mm:(10a) ', respec- 
tively. The annual runoff showed a decreasing trend, with a tendency rate of 0.67 10° m*-(10a) '. (2) The concen- 
tration degree and uneven coefficient of runoff showed a weak downward trend, and the increasing dry season run- 
off trend was evident. The seasonal distribution was more uniform, and the concentration period showed a de- 
layed trend, with a delay rate of 3.0 d-(10a) . (3) The annual runoff oscillated significantly at a scale of approxi- 
mately 44 years, and the mutation occurred in 1990. After the mutation, runoff decreased by 3.52 10° m°. The dis- 
tribution of glaciers in the basin showed a decreasing trend, whereas the vegetation cover did not show a signifi- 
cant change. (4) The contributions of climate and human activities to the runoff decrease in the Datong River 
were —17.7% and 117.7%, respectively. Precipitation was the main source of water supply in the Datong River 
Basin, and interbasin water transfer was the main driving factor for runoff reduction. 

Key words: runoff evolution; cumulative slope change rate method; climate change; human activities; Datong 


River Basin 


